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PHOTOELASTIC STRESS ANALYSIS METHODOLOGY 
AND APPLICATIONS TO NERVA STRUCTURES 


I. INTRODUCTION 


Photoelasticity is a well-recognized and continually developing technique 
for experimental stress analysis. It is unique among all experimental and 
numerical methods in regard to the complexity of the structural geometry that can 
be solved exactly. Similarly, its high resolution capability for analysis of 
rapid spatial variations of stress is unrivaled by any other experimental tech- 
nique. In spite of this technical superiority in these areas, photoelasticity 
is not routinely applied for purposes of stress analysis. It finds best applica- 
tion to specially selected problems, and the reasons for this are clear. The 
desand for Increasingly sophisticated stress analysis in recent years has been 
paralleled by developments in computer technology and numerical methods for 
stress analysis. In a great many instances the numerical "model" is of at least 
equal validity to the photoelastic model. Equally important, numerical modeling 
often lends itself better to rapid and economical production of data since 
engineering personnel customarily receive education that equips them for more 
readily absorbing the numerical rather than the experimental techniques. And, 


since numerical methods involve mathematical abstractions rather than physical 
quantities, delays and inertia in making changes in the model or its testing 
sequence are often correspondingly less. Consequently, numerical stress work 
schedules may be more reliably achieved, at least in principle. 

In view of the foregoing, it is not remarkable that stress analysis relies 
to a very great extent up^n numerical techniques. Photoelasticity retains at 
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least four important functions in industrial stress analysis, however. These 
are: (a) analysis of three-dimensional structures too complex to be modeled 

numerically within the limitations of the available programs and machines, 

(b) the experimental verification of stress predictions obtained. from new, untried 
computer programs, (c) modeling structures which are composed of two or more con- 
tacting components and in which the contact boundary conditions change with in- 
creased loading, and (d) verification of approximate methods for solving three- 
dimensional stress problems by utilizing two-dimensional numerical techniques. 
Examples where these functions have led to the employment of photoelastic tech- 
niques are plentiful in the literature. Appendix A contains some examples of 
these and other applications. Other examples are provided by References [1]* 
through [3], and Reference [12]. 

Three major applications of photoelasticity were made to NERVA structures 
components during the course of the program. A review of this experimental work 
and some of the benefits derived from it is contained in a later section of this 
report. 

II. OUTLINE OF PHQTOELASTICITY PRINCIPLES 
A. THEORY 

In this brief outline, attention will be directed to the transparent 
plastic model materials that possess isotropic elasticity and which exhibit 
temporary double refraction under stress [4], 

R = Ct (a 1 - a 2 ) (1) 


^Reference numbers are in square brackets throughout this report. 
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where R = the relative retardation of transmitted light 

C = the material stress-optic coefficient 
t = model thickness 

(o^ - Oy) = the difference of the major and minor principal stresses 


This relation is depicted graphically at the top of Figure (1) where 
a plane stressed cubical element is shown being traversed by a single ray of 
polarized light that has been resolved into in-phase components in the two planes 
of principal stress* Below, the simplest of all polariscope configurations is 
depicted, consisting of a crossed polarizer and an analyzer. The function of 
each element is described there. This polariscope provides not only the relative 
retardation caused by stress, but also direction, 0, of the principal stress with 
respect to the x-axis. 

A plane stressed cubical element oriented to the x, y coordinate 


system is subjected to the normal stress components a ^ and , and shear stress 

t . The principal stresses are related to these stress components by the static 
xy 

equilibrium requirements: 




and 0 



(4) 
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The inverse relations are: 

vj 


a 

x 


2 

a^cos 0 


c^sin 0 


(5) 


2 2 

a = a. sin 0 + a 0 cos 0 
y 1 2 


( 6 ) 


T 

xy 



sin20 


(7) 


The retardation R depicted in Figure 1 is the stress-induced spatial 
phase lag introduced in the plastic model by the difference of the principal 
stresses. Its physical origin is in the change of the refractive index accompany 
ing changes of stress in such plastics. Thus the velocity of the light ray com- 
ponent polarized in the plane is not the same as that of the other component 
beam polarized in the plane. The beams enter the element in phase but emerge 
out of phase. In terms of the engineering quantities generally employed in 
photoelastic analysis, equation (1) becomes 


n 




a 2 )t 


( 8 ) 


where n = retardation in units of multiples of the wave length of 
the monochromatic light employed, fringe order* 

f = the "s tress-fringe value”, a material constant, psi-inches 
per fringe 


*The term "birefringence” is commonly used to denote the occurrence of this 
double refraction phenomenon. Thus, n is the amount of birefringence. 
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Knowing the model thickness, the fringe order and the stress-fringe value it 
is seen that one may readily solve for the difference of principal stress at 
any point in the photoelastic model. 

Equation (8) is sufficient to determine the surface stress tangent 
to the model or slice edge, when that surface is free of contact with another 
solid body. In order to determine the true principal stresses, which are not 
necessarily coincident with the slice plane, supplementary methods must be em- 
ployed. Chief among these are the method of oblique incidence and the method 
of sub-slicing for frozen-s tressed models. [5] 

Interior stress (stress at other than boundary locations , or at 
boundary locations where contact with another solid occurs) are also investigated. 
The "shear-difference” method [6] offers the most practical approach for "frozen- 
stress" experiments, though the oblique incidence approach is also useful, in 
principle. 

B. MODEL/PROTOTYPE RELATIONS 

Stresses and deflections determined from the subscale photoelastic 
model are ultimately interpreted in terms of the full scale prototype. Model- 
prototype scaling laws emerge from dimensional analysis procedures and from con- 
sideration of the laws which govern similitude between the model and the prototype, 
[ 7 ] • These laws can be readily deduced for the simplest case of an externally 
loaded, rigid, elastic structure which behaves in a linear manner. The stress 
in a solid body which fits this description can be expressed functionally as 


o 



Z 

a 9 


z_ 

a 9 


v. 


S) 


(9) 



R. C. Sampson 
2 June 1972 


N8120R: 72-034 
Page 6 


where P = the applied force vector 

a = a characteristic dimension of the body 
x, y, z = the coordinate position of the stressed point 
v “ Poisson’s ratio, a material constant 
S = a complete definition of the shape of the body # 
If the stress is put in dimensionless form. 


P a C a’ a’ a’ V ’ S} 


( 10 ) 


the dimensionless stress is seen to be a function only of those parameters 
within the brackets. If each of these parameters is the same for the model as 
for the prototype, then their dimensionless stress at corresponding points will 
also be the same. Similar reasoning also applies to more complex conditions. 

By the same reasoning as above, the dimensionless deflections of 
either the model or the prototype are given by 


Ea6 

P 


F. (* v, S) 

6 a ’ a a* 


(ID 


Again the terms within the brackets on the right hand side must be identical 
in both the model and the prototype in order that exact predictions of prototype 
deflection may be made from observations of the model deflection. The effect 
of differences in these parameters varies considerably. Each case should be 
evaluated separately rather than relying on generalizations. 

The behavior of the viscoelastic materials employed in photo- 
elasticity becomes very complex under certain circumstances* Nearly all photo- 
elastic testing is conducted at temperatures where the physical and optical 



R. C. Sampson 
2 June 1972 


N8120R: 72-034 
Page 7 


behavior of the material is essentially independent of time. Interpretation of 
the "stress patterns" is thus greatly simplified. An outline of photoelastic 
materials behavior follows. 

C. MATERIALS BEHAVIOR 

A typical photoelastic material is characterized schematically with 
regard to birefringence (or strain) as a function of temperature and time in the 
manner illustrated in Figure (2). Briefly, there are two regions of this "contour 
map” where the photoelastic response of the material is, for all practical pur- 
poses, independent of time and/or small changes of temperature. These regions 
are those labeled ”Rubbery” and ”Glassy”. A photoelastic material of great 
practical value for plane stress modeling must have a glassy region that extends 
well above ordinary room temperature. Materials useful for three-dimensional 
"frozen-stress” analysis must have their rubbery regions occurring at tempera- 
tures low enough that surface oxidation or other thermally activated chemical 
changes do not occur to a significant degree. For most purposes, the "Transition” 
region in Figure (2) is not specifically utilized in photoelastic testing, though 
exceptions have been noted in very special investigations [8]. 

The "freezing of stress” is of special importance in photo-elasticity 
because it enables the permanent fixing of the stress-induced birefringence 
system in a three-dimensional structure. Loading fixtures are removed and the 
model dissected at leisure for examination of all the stressed locations of 
interest. In order to better describe this process, a schematic representation 
of the photoelastic material as an arrangement of springs and dashpots is helpful, 
Figure (3). This "model” is consistent with the behavior represented in Figure (2) 
when the model elements are as defined in the legend of Figure (3) . Briefly 



R. C. Sampson 
2 June 1972 


N8120R: 72-034 
Page 8 


stated. If a constant stress is applied to the material for a sufficiently long 

time, the birefringence will reach a value n where no further increases with 

max 

time will be noted, independent of temperature. However, if the material is 

first heated to its "critical 1 ' temperature, T , it will achieve the value n 

cr max 

almost instantly and maintain that value indefinitely so long as the stress and 
temperature are maintained constant. At that elevated temperature, the dashpot 
viscosities and tend toward zero and the elastic elements that predominate 
at room temperature, K^, and , no longer sustain any load. All the birefringence 
and strain is accumulated in the relatively soft "spring" K^. It must be recalled, 
however, that while Figure (3) is a functional model that behaves in the same 
general manner as the photoelastic material, the actual mechanisms involved in. 
these polymeric materials are far more complex. For that reason, functional 
models like that in Figure (3) are completely inadequate for quantitative char- 
acterization of actual materials. Nevertheless, employing this model helps to 
understand the general behavior, if not the true origins, of photoelastic testing 
principles. For example the distinction between the behavior of cross linked and 
uncrosslinked polymers is depicted in Figure (4) , the difference being in the 
values prescribed for the steady flow dashpot in Figure (3). The various 
stages of deformation, where different model elements predominate, are indicated 
by letter-designated segments on the curves of Figure (4). At time t^, the 
load was removed in this illustration, and it is seen that the uncrosslinked 
material experiences a permanent set, whereas the crosslinked material "remembers" 
its original stress free state after a period of time. 

"Stress freezing" materials are of the crosslinked type. The photo- 
elastic behavior of a typical material during the process of stress freezing is 
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shown schematically in Figure (5) . For simplicity, the behavior of a uniaxial 
tensile specimen is represented there. No essential difference of behavior is 
introduced regardless of the geometric complexity of the structural model. The 

i 

significant stages of the process are labeled on the curve in Figure (5) , with 
letters designating the terminals of each phase. It is seen that the process 
consists sequentially of a loading phase, a heating phase, a cooling phase and an 
unloading phase. The sequence of corresponding events in the functional model, 
Figure (3), can be readily appreciated. The final condition is one where spring 
element is held in an extended condition by the relatively very rigid spring 
elements and K^. Elastic energy of deformation is thus stored within the 
system, and in the actual material this storage takes place on the molecular 
scale. Thus, the cutting or destruction of a part of the model has no effect 
on those remaining elements only a few molecular crosslinks distant from the cut. 

In the manner just described schematically, the total stress system 
is permanently stored in three-dimensional structural models. Almost unlimited 
dissection of the model is permissible, since no significant release of stored 
energy occurs in the remaining material. Analysis may be initiated immediately 
or at any finite later time, and it may be as cursory or as detailed as the in- 
vestigator wishes. It is not surprising therefore that the "frozen stress" 
photoelastic technique finds far more practical application than any other method 
in photoelasticity. 

III. APPLICATIONS TO NERVA STRUCTURES PROBLEMS 

Three photoelastic investigations were conducted in the course of the 
NERVA program. These three are listed in chronological order as References [9], 
[10], and [11]. An outline of each, stating the problem and the principal test 


results, follows. 
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A. ROLLER BEARING, Reference [9] 

The roller bearing. Figure (6), was required to operate with the 
loose-fitting shaft axis misaligned as much as 0° - 30 1 w;Lth the bearing axis, 
and an unfavorable localization of load occurred on the rollers. Failures of 
the inner race were experienced in test operations. Photoelastic model tests 
were undertaken for the purpose of determining the location and magnitude of peak 
stresses. Two model tests were conducted. Test I with the cylindrical rollers 
having standard corner relief and Test II with a generous radius comer relief. 
The photoelastic plastic models were tested by the "frozen stress" technique 
in which the structurally similar model was statically loaded with appropriately 
scaled dead weights in an oven which was heated to the critical temperature, T . 
For the material employed in these tests, HYSOL 4390, the critical temperature 
was 128°C. The proper inclination of the shaft was regulated and controlled by 
adjusting the position of the weight hangars, Figure (7). 

Numerous slices were removed from the models after completion of 
the tests. Slices taken from two corresponding rollers in the two tests are 
shown in Figure (8). The stress distributions were charted at a number of 
locations on the inner race and converted to stress in the prototype by appro- 
priate similarity relations, equation (10). A comparison of the inner race 
meridional stress distribution for Test I (Figure (9)) and Test II (Figure (10)) 
is shown. 

The experimental results showed that major differences in contact 
surface stress existed in the standard and the modified roller when shaft mis- 
alignment was introduced. Modification of the bearing roller ends as in model 
Test II was concluded to offer significant improvements in the bearing stress 


concentration. 



R. C. Sampson 
2 June 1972 


N8120R: 72-034 
Page 11 


B. BUTTERFLY VALVE LIP SEAL (P/N 1117995-11), Reference [10] 

Three-dimensional photoelastic models were employed to determine 

separately the elastic stress distribution caused by seating forces and by 
pressurization forces. The lip seal gate valve prototypes provide flow control 
in a liquid hydrogen supply line. The rotating gate wipes against the edge of 
the lip-seal ring in closing, enforcing a certain radial displacement of the ring. 
When the gate valve is opened, a substantial fluid pressure differential is 
momentarily established across the seal ring. The photoelastic model test fix- 
tures were designed to model, by static testing means, these two separate loading 
conditions. 

The valve seating load test setup is shown in cross-section on 
Figure (11) and the pressure test setup in Figure (12). In each case, ’’freezing" 
of the model stresses was accomplished in the conventional manner, and meridional 
slices were removed from the seal ring for analysis. Figure (13) shows the 
photoelastic "frozen-stress" freeze pattern in a representative slice taken from 
each test. 

The analysis of these stress patterns showed that the proposed 
seating conditions will produce forces in the prototype that will induce lip 
seal ring stresses well above the yield point. Transient stresses due to the 
momentary pressure drop across the prototype seal ring, due to valve opening, 
were found to be below the yield point except very locally near the retaining 
ring clamp location. 

C. IMPELLER ROTOR, Reference [11] 

Verification of stress levels in the TPA impeller first stage rotor 
was required, due to the extreme complexity of the stress analysis problem. 

The spiral blades of the impeller posed especially severe problems in devising a 
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suitable finite element model for computation of stresses* Therefore an accurate, 
full scale photoelastic model. Figure (14), was prepared and installed in a 
spin-test fixture inside a "s tress-freezing" oven. Figure (15), Rotation of the 
model at a speed of only 1758 rpm was sufficient to induce "frozen” stresses and 
deflections that enabled a thorough analysis. Projecting the model stresses and 
deflections in terms of the prototype rotating at much higher speed was later 
accomplished through model/prototype similarity conditions similar to Equations 
(10) and (11). 

A typical "frozen-stress" slice taken from a meridional plane through 
one of the blade tips is shown at left in Figure (16). The value of the fringe 
order n is noted at several places on that photograph. These values were measured 
at all significant surface stress locations, and the boundary stress determined 
by Equation (8). In all, 41 slices were removed from various planes and the 
numerous peak surface stresses determined. Demonstrating the reversible elastic 
nature of "frozen" stresses, the slice at left in Figure (16) was reheated to 
Its test temperature and slowly cooled. Again viewing if in the polariscope, 
at right in Figure (16) , it is seen that all traces of the spin-induced elastic 
stress pattern have disappeared. 

As a follow-on to this experiment, a finite element computation model 
was developed [12] which greatly improved the accuracy of stress calculations 
above that formerly obtained on the basis of finite element models derived from 
largely intuitive representations. The improved finite element models provided 
calculated peak stresses that agreed within 10 percent, on the average, with the 
experimental results, while the maximum stresses in the models agreed within 
4 percent. Experimentation with the finite element model also indicated that 
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Poisson’s ratio has a significant effect on the stresses experienced in the 
impeller, particularly in the blade/shroud region, 
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FIG. I THE PHOTOELASTIC EFFECT 
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Figure 7 
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Figure 9 - STRESS DISTRIBUTIONS ON RACE SURFACE, TEST I 
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Figure 15, 
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1 


▼ 




Figure 16 


FROZEN STRESS" PATTERN IN SLICE NO. I (Left) 
AND AFTER ANNEALING (Right) 


R. C. Sampson N8120R: 72-034 

2 June 1972 Page 30 



2:15 p.m. Discussion 


3:00 p.m. Tea 


3:30 p.m. Polymers in Adverse Environments 

R. E. Eshenaur (Manager, Plastics Applications 
and Evaluations. General Electric Co., Kentucky, 
U.S.A. ) 

Chairman: W. F. Watson (Director, Rubber and 

Plastics Research Association of 
Great Britain, Shropshire, England) 
Panel: N. Grassie (Reader in Chemistry, 

Glasgow Univ., Scotland) 

R. S. Hagan (Manager, Plastics Lab- 
oratory', General Electric Co., Ken- 
tucky, U.S.A. ) 

J. B. Howard (Bell Telephone Lab- 
oratories, Inc., New Jersey, U.S.A.) 

L. N. Phillips ( Head of Plastics Tech- 
nology Section, Royal Aircraft Estab- 
lishment, Hants., England) 

4:00 p.m. Discussion 

5:00 p.m. Meeting Adjourns 


11:15 a.m. 
12:15 p.m. 
2:00 p.m. 


TUESDAY, JUNE 17, 1989 

9:00 a.m. High Mechanical Performance — the Engineering 
Design Problem 

R. M. Ogorkiewicz (Lecturer, Dept, of Mechan- 
ical Engineering, Imperial College, London, En- 
gland ) 

Chairman: J. H. Crate (Consultant, Plastics 

Dept., Technical Services Lab., Du 
Pont Co., Delaware, U.S.A.) 

Panel: G. M. Buehrig (Art Center College 

of Design, California, U.S.A. ) 

R. A. Horsley (Group Head, Polymer 
Physics, Shell Research Ltd., Man- 
chester, England ) 

B. W. Nelson ( Dept. Head, Plastics 
Laboratory, National Cash Register 
Co., Ohio, U.S.A.) 

D. W. Saunders ( Prof, of Polymer 
Physics and Engineering, College of 
Aeronautics, Cranfield, Bedford, En- 
gland ) 

9:30 a.m. Discussion 

10:15 a.m. Coffee 


10:45 a.m. The Problem of Polymer Selection and the Con- 
tribution of Plastics to Other Technologies 
A. II. Willboum (Director of Research, Imperial 
Chemical Industries Ltd., Plastics Div., Herts., 
England ) 

Chairman: F. R. Eirich (Dean of Research, Poly- 

technic Institute of Brooklyn, New 
York, U.S.A.) 

Panel: L. B. Allen (Manager, Materials Lab- 

oratory, International Business Ma- 
chines, Inc., New York, U.S.A.) 


2:30 p.m. 
3:15 p.m. 
3:45 p.m. 


4:15 p.m. 
5:00 p.m. 


C. E. Hollis (Technical Coordination 
Manager, BP Plastics Dept., BP 
Chemicals (UK) Ltd., London, En- 
gland ) 

II. J. Sharp (Manager, GKN Plastics 
Development Centre, London, En- 
gland ) 

G. M. Wolf ( Manager, Body Ma- 
terials Engineering Dept., Ford Mo- 
tor Co., Michigan, U.S.A.) 


Discussion 

Lunch 

Processing for High Performance 
J. F. Carley (Assoc. Prof., Dept, of Chemical En- 
gineering, Univ. of Colorado, U.S.A. ) 

Chairman: M. B. Godfrey (Managing Director, 

Thermo Plastics Ltd., Bedforshire, 
England ) 

Fanel: E. Caspar (Consultant, London, En- 

gland ) 

E. C. A. Homer (Polyolefins Tech- 
nical Service Manager, Imperial 
Chemical Industries Ltd., Plastics 
Div., Herts., England) 

J. van Leeuwen ( Head, Scientific 
Technical Dept., Plastics and Rubber 
Research Institute, TNO, Delft, Hol- 
land ) 

I. Rubin ( President, Robinson Plas- 
tics Coq)., New York, U.S.A.) 

A. Spaak ( Director, Technical Ser- 
vice, Allied Chemical Corp., New 
Jersey, U.S.A.) 

Discussion 

Tea 

New Concepts in Processing 

J. M. Goppel ( Koninklijke/Shell Plastics Labora- 
tories, Delft, Netherlands ) 

Chairman: D. V. Rosato ( Plastics World maga- 

zine, Massachusetts, U.S.A.) 

Panel: W. B. Evans (Vice President, Im- 

proved Machinery, Inc., New Hamp- 
shire, U.S.A.) 

A. Kennaway (Consultant, London, 
England ) 

G. E. Pickering (Senior Staff Asso- 
ciate, Arthur D. Little, Inc., Massa- 
chusetts, U.S.A.) 

J. Stafford (Assistant Development 
Manager, Imperial Chemical Indus- 
tries Ltd., Plastics Div., Herts., En- 
gland ) 

Discussion 
Meeting Adjourns 


SPE Journal, March 1 . 96,9 — Vol. 25 


23 


4 


8120R: 72-034 



APPENDIX A 


Structures Problem-Solving 
in the Photoelastic Laboratory 

Robert C. Sampson, Aerojet-General Corp., Sacramento, Calif. 


Photoelasticity is no longer the scientific novelty 

that formerly evoked expressions of amazement and, 
oftentimes, over-enthusiasm as to its practical applica- 
tions. While retaining its unique ability for visual 
representation of stress in structural models, photo- 
elasticity has matured and found its place due to key 
advances in materials and techniques. The employ- 
ment of photoelastic model methods at Aerojet-Gen- 
eral Corp. in support of an industrial program for en- 
suring structural integrity has been highly successful. 

In the development of space propulsion systems 
structural integrity must be maintained, but weight 
must also he reduced to a minimum. These conflicting 
requirements pose problems that can be solved in 
one of three basic ways: (a) build the structure and 
test it under actual or simulated conditions, (b) ana- 
lyze the structure by building a mathematical model 
and testing it in a computer, or (c) fabricate a true 
scale physical model of the structure and measure its 
structural performance in closely-controlled labora- 
tory conditions. In practice, methods (a) and (h) 
provide the major share of the required solutions, but 
the third alternative can be extremely useful in cer- 
tain special cases. It is this category of special prob- 
lems that are described here. 

Six examples have been selected from test records 
to illustrate the variety and scope of the problems 


that have found most effective treatment by the 
model techniques. While these examples necessarily 
reflect a preoccupation with the problems of the 
rocket motor industry, a similar variety of applications 
will be found in any industry where strength and 
weight requirements are in severe conflict, and where 
there are heavy economic penalties for structural 
malfunction. 

The selected examples illustrate the general areas of 
utilization: (a) analysis of a complex structure, (b) 
computer methods verification, (c) materials research, 
(d) structural design support, and (e) two examples 
of special technique development to match current 
problems. 

Analysis of a Complex Structure 

Problem: A pressurized rocket motor chamber is too 
complex in shape to analyze the stresses by mathe- 
matical methods. How can these stresses be deter- 
mined? 

Solution: Make a true scale model from photoelastic 
plastics, then “freeze” the stresses by heating and cool- 
ing the model while the internal pressure is held con- 
stant. Cut thin slices from the model and analyze the 
frozen stresses in a polariscope. 

Example: The intersecting cylinders of the model 
( Figure 1 ) were cast integral to finish dimensions 
from an epoxy resin mixture in an acrylic plastic mold. 
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Figure 1 . Rocket chamber photoelastic model. 



A unique technique is described for visual representation of stresses in structural models 


Stresses were frozen under internal fluid pressure 
( Test 1 , Figure 2) and fluid pressure combined with 
a mechanically exerted thrust acting on the flanged 
top and reacted at the integral pad on the bottom 
(Test 2 , Figure 2). The transverse section frozen 
stress patterns (Figure 2) gave all the information 
needed to accurately depict the surface stress distri- 
butions in that plane. The numerical values of fringe 
order, n, are noted at several locations on Figure 2. 



Figure 2. Frozen stress patterns in transverse sections. 


These values are directly proportional to the stress at 
each point. Stresses were also determined on a number 
of other sections, including the line of intersection of 
the two cylindrical surfaces (Figure 3). 

Comment: A stress problem too difficult for theo- 
retical analysis was solved bv a combination of good 
model-making technique and the frozen-stress photo- 
elastic method of stress analysis. 
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Figure 4. Photoelastic (right) vs 
computer (left) solution of thermal 
shrinkage problem in case-bonded 
grain. 


Verifying a Computer Analysis 

Problem: Computer solutions of a class of axisym- 
motric loading problems in case-bonded solid rocket 
motors had been obtained. It was essential to demon- 
strate the accuracy of the solution before it was rou- 
tinely employed in evaluating prototype motor struc- 
tural integrity. 

Solution: Prepare a three-dimensional photoelastic 
model and freeze the stresses under carefully con- 
trolled conditions. Analyze the stress in slices taken 
from the model and compare the results with the 
computer stress predictions. 

Example: The problem of thermal shrinkage in 
short, axisymmetric, case-bonded rocket motors had 
recently provided a numerical solution. Using the ex- 
perimental technique of the problem on thermal 
shrinkage in case-bonded propellant grains, a photo- 
elastic solution to this was obtained. The computer 
data were used to construct a map of half-order bire- 
fringence in a meridional plane of the short, case- 
bonded grain with a concentric circular bore. This 
artificial stress pattern is shown compared to the real 
photoelastic pattern in Figure 4. It is seen that very 
good over-all agreement of theory and experiment is 
indicated. A more explicit comparison was seen in the 


individual principle stresses along the central radial 
section, giving excellent quantitative accuracy. 

Comment: Confidence in the inherent accuracy of 
the numerical calculation method was firmly estab- 
lished. 

Materials Research Support 

Problem: Composite solid propellants consist essen- 
tially of rigid particles of oxidizer embedded in a rub- 
bery fuel matrix. Mechanical failure begins at the 
microstructure level in the binder matrix between 
rigid particles. Investigate the mechanical aspects of 
failure initiation between embedded rigid particles. 

Solution: Construct superscale photoelastic models 
of plane sections having rigid disc particles embedded 
in a soft matrix. Analyze the stresses between adjoin- 
ing particles and predict the probable failure mechan- 
ism. Test these predictions with experiments on arrays 
of spherical glass beads in a rubbery binder. 

Example: Rigid epoxy discs in a matrix of soft 
photoelastic rubber exhibited stress patterns like that 
in Figure 5 (left) when a uniaxial pull was exerted. 
The nearest neighboring discs, aligned with the load, 
exert a severe stress concentrated on the binder, as 
shown by the curves of stress distribution. Figure 5 




Figure 5. Rigid discs embedded in 
rubber— photoelastic pattern (left) 
and stress distributions (right). 
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Frame No. 

0 

1 

5 


15 

17 

21 


Figure 6. Progressive damage in a composite model. 


(right). Failure may equally likely initiate at the 
disc/binder interface bond or in the interior of the 
binder ligament. These conclusions were tested by 
observing the failure behavior of 5 x 9 arrays of glass 
balls embedded in a soft rubber binder. Sequence 
photographs from one such test ( Figure 6) show 



Figure 7. Photoelastic test specimen of a propellant 
grain. 


early initiation of cohesive failures between neighbor- 
ing beads aligned with the load. Frame 1. Marked 
changes of the stress-strain characteristic were noted 
near Frames 5 and 17, and a local weakness leading 
to eventual failure is evident in Frame 21. Similar 
experiments with different binders showed how the 
structural performance of the composite is most 
strongly affected by the binder properties, and that 
different failure processes result from different binder 
properties. 

Comment: Superscale model techniques are often 
very useful in delineating the mechanical principles of 
behavior in microscopic composites. Insights are 
gained that might be obscured in ‘Teal” composites. 

Stress Data in Support of Structures Design 

Problem: Grain design engineers must often weigh 
the effects of grain geometry changes on both the 
stress and combustion efficiency of the design. Design 
data are needed that relates the stress concentration 
at the “star points” of a grain cross section to its geo- 
metric parameters. 

Solution: Design, build and test a series of two-di- 
mensional photoelastic models that encompass a wide 
range of the important geometric parameters. Expand 
the useful scope of the test data by cross-plotting and 
interpolating. 

Example: A typical 2-D model is shown in Figure 7. 
The frozen stress pattern was induced by the thermal 
shrinkage stress-freezing process described in the 
problem of thermal shrinkage in case-bonded propel- 
lant grains. Stress concentration data (maximum fillet 
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stress per unit external pressure) obtained from a 
series of models of this type are shown in Figure 8. 
These curves allow the designer to assess the conse- 
quences of a wide variety of geometric changes, with 
good engineering accuracy. 

Comment: In plane structures with few shape 
parameters, extensive design data can often be ob- 
tained at moderate cost. 


Propellant Strain Measurement Technique 
Development 

Problem: The behavior of solid propellants in the 
presence of high strain concentrations is not readily 
predicted. Develop a convenient means for measuring 
surface strain distributions on actual propellant speci- 
mens. 

Solution: The moir& grid-interference method was 
adapted to this purpose. Investigation showed that 
the most effective approach utilized a grid pattern 
printed on the surface of the test specimen. Grid 
photographs of the loaded specimen, obtained at 1:1 
magnification, contain all the information needed to 
determine the surface strain distribution. A novel 
“grid-shift” technique was developed to speed the 
accurate analysis of the moire fringe data. 

Example: A slotted, wedge-shaped propellant test 
specimen closely represents the stress concentration 
in a solid rocket grain. A specimen of this kind with 
printed surface grids of 100 lines/in. was photo- 
graphed, and quantitative analysis of the moire grid 
photographs ( Figure 9) showed that the strain con- 
centration factor increased repidly as failure at the 
notch was approached. 

Comment: Moire methods allow detailed strain 
analysis on real solid propellant materials tested in 
realistic environmental conditions. 



Figure 8. Effect of slot-end configuration on maximum 
fillet stress. 


Thermal Shrinkage in Case-Bonded Propellant 
Grains 

Problem: Solid propellants expand and contract 
with temperature approximately ten times more than 
the rigid outer case materials, and severe interaction 
between the case and the grain is encountered. How 
can these stresses and strains be evaluated lor complex 
three-dimensional motors? Conventional stress-freez- 
ing photoelastic methods are not suitable for this class 
of problem. 



Figure 9. Moire displacement fringe pattern in propel- 
lant test sample. 



Figure 10. Casting mold and finned cylindrical core for 
3-D grain model. 
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LONGITUDINAL SLICES 


Figure 11. Frozen stress slices from 3-D grain model thermal shrinkage test. 


Solution: A novel stress-freezing technique was de- 
veloped to permit photoelastic simulation of propellant 
shrinkage conditions. In this method, a true-scale 
model is cast to shape in liquid epoxy resin and cured 
at gradually elevated temperature while bonded to 
the subscale model case. Cooling the composite model 
structure induces the same stress condition experi- 
enced by the motor. Differential shrinkage stress is 
frozen in and is accurately resolved in a quantitative 
manner, using special model materials and calibration 
techniques. 

Example: The subscale metal case and cores are 
shown in Figure 10. The large fins attached at the 
forward end of the cylindrical aluminum core are 
made from Wood’s metal to facilitate easy removal 
after solidification of the cast liquid resin. After the 
stress-freezing cycle is completed, the metal case is 
slit and peeled away. Thin slices taken at various 
sections and viewed in polarized light ( Figure 1 1 ) 
give all the information needed for detailed stress 
determinations. 

Comment: Differential contraction in verv complex 
composite structures was made amenable to experi- 
mental model analysis by devising a new, three-di- 
mensional stress-freezing technique. This type of 
problem has widespread significance in solid propel- 
lant rocket motor structural integrity. 

Conclusion 

These representative examples illustrate that the 
photoelastic laboratory is far removed from being a 
routine testing facility. On the contrary, the ability to 
modify existing techniques or devise new ones to cope 
with non-routine problems is readily seen to be its 
most valued asset. Maintaining and improving this 
capability is one of the laboratory’s primary functions. 
Consequently, the laboratory maintains a steady re- 
search effort aimed at improving model materials, ex- 


perimental techniques and instrumentation. Only in 
that way can a suitable repertoire of techniques be 
made ready for the new problems that periodically 
appear in an advancing industry. 
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